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Abstract. The disruption of stellar systems, such as open clus- 
ters or stellar complexes, stands out as one of the most rea- 
sonable physical processes accounting for the young moving 
groups observed in the solar neighbourhood. In the present 
study we analyse some of the mechanisms that are important 
in the kinematic evolution of a group of unbound stars, such 
as the focusing phenomenon and its ability to recover the ob- 
served moving group's velocity dispersions, and the efficiency 
of disc heating and galactic differential rotation in disrupting 
unbound stellar systems. Our main tools used to perform this 
analysis are both the epicycle theory and the integration of the 
equations of motion using a realistic gravitational potential of 
the Galaxy. 

The study of the trajectories followed by stars in each of 
the Pleiad es moving group substructures found by Asiain et 
al. (1999) allows us to determine their stellar spatial and ve- 
locity distribution evolution. The kinematic properties of these 
substructures are compared to those of a simulated stellar com- 
plex which has evolved under the influence of the galactic 
gravitational potential and the disc heating. We conclude that 
a constant diffusion coefficient compatible with the observa- 
tional heating law is able to explain the velocity and spatial 
dispersions of the Pleiades moving group substructures that are 
younger than ~ 1.5 • 10^ yr. 

Key words: Stars: kinematics - Galaxy: kinematics and dy- 
namics - { Galaxy:) solar neighbourhood 



1. Introduction 

A "supercluster" can be defined as a group of stars gravitation- 
ally unbound that share the same kinematics and may occupy 
extended regions of the Galaxy. A "moving group" (MG here- 
after) or "moving cluster", however, is the part of the super- 
cluster that can be observed from the Earth (see, for instance, 
Eggen 1994). 

Independent of the physical process that leads to the for- 
mation of superclusters, there are basically two factors act- 
ing against their persistence in the general stellar background. 



First, when a group of unbound stars concentrated in the phase 
space has a certain velocity dispersion, the galactic differential 
rotation tends to spread them out very quickly in the direction 
of the galactic rotation (e.g., Woolley 1960). Second, observa- 
tions show that the velocity dispersion of disc stars is increased 
with their age, which is usually interpreted as the result of a 
continuous process of gravitational acceleration that may be 
produced by different agents (e.g., Lacey 1991). This second 
factor, generally referred to as disc heating, also disperses stars 
very efficiently. It is therefore striking to verify that some of the 
classical MGs are several 10^ yr old. Some authors have over- 
come part of this problem by assuming that the velocity dis- 
persion of these stellar gro ups m ust be very small - e.g. Eggen 
( |1989| ), Yuan & Waxman ( ll977| ), So derblom & Mayor ( |1993| ). 
Nevertheless, some recent studies (Chen et al., 1997; Sabas, 



1997; Figueras et al., 1997; Chereul et al., 1998 



Asiain et al.. 



1999) do not support such a small dispersions. 

There have been several attempts to explain the origin of 
superclusters. Perhaps one of the most widely accepted expla- 
nation is the "evaporation" of the outermost stellar component 
(or corona) of open clusters (e.g. Efremov 1988). Over time, 
open clusters are disrupted by the gravitational interaction with 
massive objects in the Galaxy (such as giant molecular clouds), 
and as a result the open cluster members fill a long tube in 
space. The part of this tube that we can observe should have 
a small velocity dispersion - e.g. Weidemann et al. (1992), 
Eggen ( 1994|). Alt hough based in old stellar evolutionary mod- 
els by Iben ( |1967[ ), Wielen ( |1971| ) found that about 50% of the 
open clusters disintegrate in less than 210* years, a result sup- 



ported by posterior studies (lerlevich, 1987; Iheuns, 1992) 



If superclusters come from the evaporation of open clusters' s 



coronae, then Wielen's (1971) results indicate that most super- 
clusters should be young or the last tracers of former clusters. 

MGs may also be produced by the dissolution of larger 
stellar agglomerations, such as stellar complexes or fragments 



of old spiral arm 


s (IWoolley, 1960 




Wielen, 1971 




Yuan, 1977 


|31megreen, 1985 


; Efremov, 1988; 


Comeron et al. 


1997 


). Thus, 
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a MG could be a mixture of stars coming from different open 
clusters' coronae and disrupted associations that share the same 
origin (and motion). This process of MG formation includes 
the first hy pothesis mentioned here. Alternatively, Casertano et 
al. ( 1993 ) proposed that MGs are open clusters that evolved 
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under the gravitational influence of a large local mass that sur- 
rounds and traps them. However, it seems difficult to spec- 
ify the nature of such a large mass. Finally, it has been pro- 
posed that superclusters could actually be made of stars trapped 
around stable periodic orbits (e.g. Miillari et al. 1994; Raboud 
etal. 1998). 

In this paper we assume that superclusters occupied a small 
volume in the phase space in their first stages, when they be- 
came gravitationally unbound. From that point on, they evolved 
under the influence of the gravitational potential of the Galaxy, 
as it would be the case in the classical picture of stellar evap- 
oration from open clusters' coronae. We then analyse different 
aspects related to the evolution of MGs. 

First, in Sect. ^ we use epicycle approximation to study the 
ability of the "focusing phenomenon" (Yuan, 1977) to period- 



ically group MG stars when no disc heating is considered. In 
Sect, ^we show how to determine the stellar trajectories from 
an analytic expression of the galactic potential. The disc heat- 
ing effect on moving groups' properties can be considered by 
introducing random perturbations to the velocity components 
of stars. These trajectories allow us to study the evolution of 
unbound stellar systems independent of the local galactic prop- 
erties (in contrast with epicycle approximation). 

We focus our study on the origin and evolution of the 
Pleiades moving group, since it is the youngest moving group 
in the solar neighbourhood and, therefore, we may find it easier 
to recover its past properties with certain confidence (Sect. 
We consider the different Pleiades substructures found in Asi- 
ain et al. (1999, Paper I), and estimate their kinematic age and 
galactic position at birth. Finally, in Sect. || we simulate a stel- 
lar complex and determine the trajectory of its stars to study its 
disruption over time. The disc heating effect on the evolution 
of the stellar complex is evaluated. 



2. The focusing phenomenon 

The epicycle approximation allows us to study, in an analyti- 
cal way, the evolution of an unbound system of stars under the 
influence of the galactic potential. Under such approximation, 
the equations of motion of stars can be expressed as: 



v' 

c 



Ca COs(j/t + -0) 



^sin(.i- 



(1) 



in the coordinate system (i^', 77', (') centered at the current po- 
sition of the Sun. points towards the galactic center (GC), 
r]' is a linear coordinate measured along a circumference of ra- 
dius Ro (galactocentric distance of the Sun) and positive in the 
sense of the galactic rotation (GR), and points towards the 
north galactic pole (NGP). k is the epicycUc frequency, v is the 
vertical frequency, t is the time (=0 at present) and luq is the an- 
gular velocity of the Galaxy at the current position of the Sun. 
£,a,£,'b'VajCa''t' ™d tp are integration constants related to the 



current position and velocity of a star {^'^, 77^, f]'^, Q'^) as: 

2B 
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Sa 
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^° ^ 1/2 



— arctan 



arctan 



(2) 



where the B and A are the Oort's constants. In|A| we use Eqs. [l| 
to evaluate the evolution of space and velocity dispersions of a 
group of stars with low peculiar velocity with respect to their 



Regional Standard of Rest (RSR). Eqs. A.. 2 show that disper- 



sions in position and velocity components oscillate around con- 
stant values, except for the azimuthal coordinate 77', whose dis- 
persion (Jn' increases with t. After a few 10^ yr, this increase is 
dominated by the secular terms for typical position and velocity 
dispersions of young stellar groups, and it can be approximated 
by the linear relationship 
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1/2 



t 



(3) 



If we consider a sample of stars with cr^' w pc and 
(T,j' « km s^^, then all dispersions in position and veloc- 
ity would oscillate around their mean values with an epicycle 

27r 

frequency n. Stars would meet every At — — (at the position 
of the Sun At w 1.5 • 10* yr), a fact that has been referred to 



as "focusing phenomenon" (Yuan, 1977). However, when we 
apply Eq. || to a more realistic case where dispersions are close 
to those of the open clusters and associations, i.e. a few parsecs 
in <Tj/ and ^1-2 km s^^ in cr^/ , we obtain cr^/ k, 1 kpc after 
atsa5-1010*yr-a period of time shorter than the age of 
some moving groups detected in the solar neighbourhood (e.g. 
Chen et al. 1997; Paper I). Thus, galactic differential rotation 
disrupts very efficiently unbound systems of stars. Nonetheless, 
if the distribution of stars in -q' is supposed to be gaussian after 
the disruption of the stellar system, an important proportion of 
the original sample will be still concentrated in space after a 
long period of time; for instance, for cr^' = 1 kpc we can find 
still ~ 24 % of the initial stars in a region 600 pc long in 77^ 
yr 



3. Stellar trajectories 

The epicycle approximation holds for stars that follow almost 
circular orbits around the center of the Galaxy. Although valid 
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for most of our stars, some of them, with high peculiar veloci- 
ties with respect to the Local Standard of Rest (LSR), perform 
radial galactic excursions of a few kpc in length. The galac- 
tic gravitational potential changes significatively at different 
points of their trajectories, and so the first order approximation 
is no longer valid. In addition, the vertical motion is poorly de- 
scribed by a harmonic oscillation as stars gain certain height 
over the galactic plane. For a more rigorous analysis, we need 
to use a realistic model of the galactic gravitational potential. 
Stellar orbits will be determin ed fr om this model by integrating 
the equations of motion (Sect. 'iA). Moreover, the addition of a 
constant scattering in this process will allow us to a ccount for 
the disc heating effect on the stellar trajectories (3.2). 



3.1. The galactic potential 

Expressed in a cartesian coordinate system (^, 77, Q centered 
at the position of the Sun|^ and rotating at a constant angular 
velocity loq. 



\ 




(4) 



the equations of motion of a star, assuming that the gravita- 
tional potential of the Galaxy ^q{R, 9, z; t) (in cylindrical co- 
ordinates centered on the GC) is known, are: 



9£g 
di] 



ujI{Ro - ~ 2wo?7 



(5) 



A fourth order Runge-Kutta integrator allows us to numerically 
solve these equations when $0 is known, obtaining the trajec- 
tory of the star To get a realistic estimation of the galactic grav- 
itational potential we decompose it into three parts: the general 
axisymmetric potential $aS; the spiral arm $sp, and the cen- 
tral bar $B perturbations to the first contribution, i.e. 

*G = $AS + *Sp + $B ■ 



We adopt the model developed by Allen & Santillan (1991) 
for the axisymmetric part of the potential, ^as{R, z), because 
of both its mathematical simplicity - which allows us to deter- 
mine orbits with a very low CPU consumption - and its updated 
parameters. The model consists of a spherical central bulge and 
a disk, both of the Miyamoto-Nagai ( 1975 ) form, plus a mas- 
sive spherical halo. This model is symmetrical with respect to 
an axis and a plane. The authors adopted the recommendations 
of the lAU ( Kerr & Lynden-Bell, 1986 ) for the galactocentric 
distance of the Sun (Ro = 8.5 kpc) and the circular velocity at 
the position of the Sun {Qo= 220 km s^^). The Oort constants 



^ points towards the GC, rj towards the sense of the GR, and 
towards the NGP 



derived from this model are well within the currently accepted 
values. 

Spiral arm perturbation to the potential is taken from Lin 
and associates' theory (e.g. Lin 1971 and references therein), 
that is: 



$sp(R, 0; t) = Acos{m{npt -e) + (P{R)) 



(6) 



where 

A = 



(i?o'i^o)^/rotani 



m 

m f R 
In -— 

tan i \Ro 



A is the amplitude of the potential, fro is the ratio between 
the radial component of the force due to the spiral arms and 
that due to the general galactic field, flp is the constant angular 
velocity of the spiral pattern, m is the number of arms, i is 
the pitch angle, (j> is the radial phase of the wave and 0o is 
a constant that fixes the posit ion of the minimum of the spiral 
potential. According to Yuan ( 1969 ), we adopt fro — 0.05 (this 
value has been confirmed in a recent study by Fernandez, 1998) 
and flp = 13.5 km s^^ kpc^^. 

W e consider a classical two arm pattern (to — 2) (Yuan, 
1969; Vallee, 1995 ). If we assume that the Sagittarius arm is lo- 
cated at a galactocentric distance Rsag = 7.0 kpc, and that the 
interarm distance in the Sun position is AR =3.5 kpc as sug- 
gested by observations on spiral arm tracers (Becker & Fenkart, 
1970; [Georgelin & GeorgeUn, 1976| ; [Liszt, 1985| ; Kurtz et al., 
1994), then the pitch angle can be determined to be: 



tan i = — ln(l + Rg^^ AR) . (7) 

7? 35 (close to the Yuan's (1969) 



27r 



For TO = 2, we obtain i 
value i = 6?2). 

For t — and = 0°, the adopted value (/)o leads to a mini- 
mum in the potential at the observed position of the Sagittarius 
arm (R = Rsag — 7.0 kpc, I = 0°). (po can thus be expressed 
as: 



.In 

Ro 



tanz 



(8) 



For i = 7?35 and to = 2 we obtain 00 = 0.131 rad. 

The central bar potential we use here is, for simplicity, a tri- 
axial ellipsoid with parameters taken from Palous et al. (1993) 
i.e. 



GMba 



Ibar 



"'bar y2 
^bar 



^bar 1 
'^bar 



1/2 ' 



(9) 



where x = Ro c os(6>— f^gt— 6*0) and j/ = Rp sin (6>— Slgt— 6*0), 
with 6*0 =45° ( fWhitelock & Catchpole, 1992| ). abar, har and 
Cbar are the three semi-axes of the bar, with qbar its scale 

^bar 1 ^bar 



length, and with 



1 



bbar 0.42 

Cbar 0.33 
5 kpc. The adopted total mass of the bar, M^ar^ 
and its angular velocity, ilg, is 70 km s^^ kpc~^ (Binney et 



and qbar = 
is 10^ M0, 
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al., 1991). Although most of the parameters that define the bar 
are very uncertain, the effect of the bar on the stellar trajec- 
tories becomes important only after several galactic rotations, 
which requires a length of time greater than the age of the stars 
considered in our study. 

3.2. Disc heating in the stellar trajectories 

The observational increase in the total stellar velocity disper- 
sion ((t) with time (i), or disc heating, can be approximated by 
an equation of the form: 



(10) 



where cTo is the dis persion at birth and Cy the "apparent dif- 
fusion coefficient" (Wielen, 1977). The constants n, a a and 



Ci, give one important informatio n on the phys ical mechanism 
responsible for the disc heating ( |Lacey, 1991 ; Fridman et al., 
1994). 

Our sample of B and A main sequence type stars, de- 
scribed in Paper I, allows us to determine an accurate and de- 
tailed disc heating law for the last 10^ yr, given the quality 
and uniformity of our ages and the size of this sample (2 061 
stars). A standard nonlinear least-squares method (Levenberg- 
Marquardt method) has been applied to fit the heating co- 
efficients to our data (Fig. |l]). In this way we obtain Uo ~ 
12 km s^^, n « 5 and C„ « 0.01 (km s^^)"yr^^, similar 
to Lacey's (1991) results for n — b, i.e. CTo = 15 km s^^ and 
Cv — 0.01 (km s~^)^yr~^. When we fix n to 2, we obtain 
Co ~ 15 km s^^and « 5.6 • 10^^ (km s^^)^yr^^, which 



is again in agreement with Lacey's ( |1991| ) (cto = 15 km s ^ 
and a, = 5 • lO"'^ (km s-i)2yr-'^Tand Wielen's ( [l977| ) 
((To = 10 km s-i and C„ 6 • 10"^ (km s-i)2yr-Tfits 
with n = 2. 

Our accurate observational heating law (circles in Fig. [l|) is 
far from being smooth. Two special features can be observed 
on this law: first, the velocity dispersion shows a steep increase 
during the first ~ 4 • 10^ yr, probably due to the phase mix- 
ing of young stars. After this point, this increase becomes less 
pronounced. Second, an almost periodic oscillation seems to be 
superimposed over the "continuous" heating law. This oscilla- 
tion (period « 3 • 10^ yr) could be the signature of an episodic 
event (see for instance Binney & Lacey 1988, Sellwood 1999). 

Since we do not know the details about the mechanisms 
that perturb the stellar orbits and cause the disc heating ef- 
fect, we assume for simplicity that the accelerating processes 
acting on a given star can be approximated by a sequence of 
independent and random perturbations of short duration (Wie- 
len, 1977). In order to evaluate the magnitude of these per- 
turbations we generated a sample of 500 stars located around 
the position of the Sun at i = 0, and uniformly distributed 
in a 100 pc/side cube. Initially, these stars are moving with 
the same velocity as the LSR, with an isotropic dispersion 
in velocity ao- During the process of orbit determination, the 
velocity of the star is instantaneously changed by /S.v every 
A yr This At? is normally distributed around zero with a (con- 
stant) isotropic dispersion <7h in each component. The equa- 




Fig. 1. Fit of Equation |l^ to our data (see text), using stars 
closer than 300 pc from the Sun, with 100 stars per bin. Solid 
line: CTo, and n are fitted; dashed line: Go and Cv are fitted, 
whereas n is fixed to 2. The error bars on points include only 
statistical uncertainties, calculated as ta — o'/V2]V, where N 
is the number of stars per bin 



30 
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Fig. 2. Velocity dispersion vs. time, determined by integrat- 
ing the equations of motion of a simulated sample when an 
isotropic diffusion is applied at each time-step. Dashed line: 
observational fit to disc heating with a fixed slope n=2 



tions of motion are integrated using only the axisymmetric part 
of the galactic potential ('I' as)- This allows us to introduce the 
heating effect only through the parameter cr/j, avoiding possi- 
ble redundances, such as the scattering of disc stars by spiral 
arms. The best approximation to the observational heating law 
is obtained when = 1.45 km s~^ for At = 10^ yr, and 
(To ~ 15 km s^^. It can be demonstrated that Gh and At are re- 
lated with the "true diffusion coefficient" D introduced by Wie- 

2 

— ^. On the other hand C,, 
At 



len ( |1977| )just as I> = On the other hand C„ = 2.95 D, 
as demonstrated by Wielen ( 1977 ) using the epicycle theory. 
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We obtain in this way an independent estimation of C^, i e. 
C„ = 2.95 — ^ = 6.2 • 10"'^ (km s-i)2yr-i. This value 



10^ 



is in an excellent agreement with our previous estimation. In 
Figure ^ the evolution of the total velocity dispersion for this 
simulated sample is compared with the fit of Eq. |l^ to the ob- 
servational data when n = 2, showing an excellent match. This 
procedure will be used in Sect. H to simulate the evolution of a 
stellar system under the influence of the disc heating. 

4. The Pleiades Moving Groups 

In Paper I we used a method based on non-parametric den- 
sity estimators to detect MGs among a sample of 2061 B 
and A main sequence type stars in the four-dimensional space 
(U,V,W,log (age) ). We used HIPPARCOS data as well as radial 
velocities from several sources (see Paper I for the details) to 
determine the stellar spatial velocities, and Stromgren photom- 
etry for the stellar ages. Tables 1 and 2 in Paper I show the main 
properties of the MGs when separated in the (U,V,W,log (age)) 
and (U,V,log (age)) spaces respectively. Since the W-velocity 
component is less discriminant than the other three variables, 
we will focus our study on the MGs in Table 2 of Paper I (Ta- 
ble |l| here). In particular, as akeady mentioned in the Sect. |l], 
we will deal with those MGS whose velocity components re- 
semble those of Pleiades open cluster. 

Table 1. Number of members {N), velocity components and 
ages (along with their dispersions) of the Pleiades MGs found 
in Paper I 



Moving 


M 


U 


V 


W 


T 


Group 




[kms^i] 


[kms-1] 


[kms-^] 


[10** yr] 


Bl 


34 


-4.5(4.7) 


-20.1(3.3) 


-5.5(1.9) 


0-2(0.1) 


B2 


75 


-10.7(5.3) 


-18.8(3.7) 


-5.6(2.2) 


0.6(0.2) 


B3 


50 


-16.8(5.1) 


-21.7(2.7) 


-5.6(4.6) 


3.0(1.2) 


B4 


53 


-8.7(4.8) 


-26.4(3.3) 


-8.5(4.7) 


1-5(0.5) 



It is interesting to com pare th ese results with those recently 
obtained by Chereul et al. ( 199S| ), who also found similar sub- 
structures in the Pleiades and other MGs by means of a wavelet 
analysis performed at different scales. The velocity dispersions 
of their substructures are quite a bit larger than those found for 
classical MGs, in agreement with our results. In their analysis, 
they did not consider the stellar age as a discriminant variable, 
which prevents them from detecting those substructures that 
are strongly defined in age but not so well defined in the ve- 
locity components. Another important difference between both 



methods is that Chereul et al. (1998) did not determine pho 



tometric ages for AO to A3 stars, since no reliable metallicity 
is available for them. Instead, they computed a "paliative" age 
that produces an artificial peak in age of 10^ yr, and a lack of 
other you ng star s (up to ^ 5 • 10'* yr). As mentioned in Asi- 
ain et al. ( 1997 ), a metallicity Z = 0.02 is representative, in 
a statistical sense, of (normal) A type stars. Using this value. 



we did not observe any lack of stars in the young part of the 
age distribution (Paper I). On the other hand, since we do not 
use F type stars in our study because of the high uncertainties 
involved in the process to determine their ages, our data do not 
allow us to confirm the existence of the ^ 10^ yr old Pleiades 



substructure found by Chereul et al. (1998 ). 

Using the numerical integration procedure described in 
Sect. [3.1| , and the mean properties (nuclei) of the Pleiades sub- 
structures found in Paper I (Table [l|), we have computed the 
trajectories of these substructures from the present up to the 
moment they were born (Fig. |]). This latter age is defined as 
the average age of the MGs constituent members. The youngest 
group, i.e. Bl, is composed of Scorpio-Centaurus (Sco-Cen) 
OB association members (Paper I), and it was born in the inter- 
arm region. In Sect. \A we study the evolution of this group. 
Since B 1 is still too young to be affected by the disc heating 
effect or phase space mixing, we can determine its kinematic 
age with som e co nfidence. The B2 group is considered sepa- 
rately in Sect. \.2 This is also quite a young group and contains 
accurate information on some of the closest associations. The 
birthplace of the older groups, i.e. B3 and B4, is close to a min- 
imum of the spiral arm potential, which seems consistent with 
their being born around this structure. Details on their spatial 



and velocity evolution are given in Section 4.3. 

In recent studies based on the velocity fiel d of C epheids, 
Mishurov et al. ( |1997[ ) and Mishurov & Zenina ( |1999| ) obtained 
a set of spiral arm parameters which clearly differ from those 
adopted here (e.g., in Mishurov & Zenina ( |1999| ), 00 = 142° 
and Q,p — LOo ~ 0.5 km s^^ kpc^^). When considering these 
new parameters, all Pleiades MG substructures turn out to be 
also placed, at birth, around the spiral arms. In this particular 
case, the Pleiades substructures appear, at birth, concentrated 
in a small area of the Galaxy. We also tested that the main con- 
clusions of the current paper does not critically change with the 
selection of any of this two spiral arm patterns, as expected. 

4.1. Scorpio-Centaurus association 

Because of the short age of these stars and the quality of our 
data we are able to determine quite precisely these stars' kine- 
matic age, defined as the time at which they were most con- 
centrated in space - assuming that they are gravitationally un- 
bound. We consider here only those stars in Bl that are con- 
centrated in space (see Fig. 7 in Paper I). With the exception of 
HIP84970 and HIP74449, aU of these stars were classified as 
members of Sco-Cen association by de Zeeuw et al. (de Zeeuw 
et al., 1999) ^ Their position in the galactic and meridional 
planes as a function of time are shown in Figure ^ It is quite 
evident that between around 4 and 12 Myr ago these stars were 
closer to each other than they are at present. However, obser- 
vational errors produce additional dispersions on velocity and 



- Though Bl is considered as a single group because of the few 
members it contains, most of its stars belong to three neighbouring as- 
sociations: Upper Scorpius (U S), Upper Centaurus Lupus (UCL), and 
Lower Centaurus Crux (LCC) (|de Geus et al., 1989|; pTaauw, 199 1|; d e 
Zeeuw et al., 1999). 
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Fig. 3. Trajectories of the MG's nuclei in Table [ijbackwards in 
time, from present (t = 0) until their mean age. The dashed 
lines represent spiral arms, as defined in Sect. The refer- 
ence system is rotating with the same angular velocity as the 
spiral arms (Hp) 



position, and therefore the maximum spatial concentration we 
find is shifted to the present. Figure | also shows us that the 
last intersection of the Sco-Cen association with the galactic 
plane was between 8 and 16 million years ago. 

The evolution of the errors or dispersions in the stellar po- 
sition and velocity can be calculate d mo re precisely by means 
of the epicycle approximation (Eqs. ). We now consider the 
fact that the observed dispersion in position inside a given MG 
(cp.o&s, where p is the distance from the LSR to a given star) 
can be decomposed into two parts, i.e. 



(T 



p,oos 



(7 



p.int 



a 



p,err 



it). 



(11) 



where 

<^p.int{t) is the MG intrinsic dispersion at the time t, and 
Cp^errit) is the MG mean observational error 

Since both a-p^obs and Cp^err can be easily determined at dif- 
ferent epochs (crp obs(i) can be calculated by integrating the 
stellar orbits until time t, and ap^crr (t) can be propagated from 



present using Eqs. |A.2| ), we can also derive the time tmin at 
which ap^int (t) was minimum. In Figure ^ we show the evolu- 



tion of cr 



p,obs 1 ^ p.crr 



and a 



p.int 



We observe a clear minimum in 



o'p,mt around tmin ~ 9 • 10 yr, which corresponds to Bl kine- 
matic age. At the minimum spatial concentration we find the 
intrinsic spatial dispersions are a^' = 20 pc, ct^' = 22 pc, and 
(T(^' = 15 pc, smaller than their current values of a^' = 30 pc, 
cr^' = 33 pc, and a^^' = 16 pc. Also, their intrinsic velocity dis- 
persions were smaller at tmin, and < 2 km s^^ in each compo- 
nent. 

There is a clear discrepancy between the kinematic and the 
"photometric" age of B 1 . This could be due to several reasons. 
First of all, most of the stars in Bl are massive, and their at- 
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Fig. 4. Spatial distribution of the 27 stars in B 1 sub-group that 
are spatially concentrated in space as a function of time. Num- 
bers in the upper right corners indicate the time in Myr 



mospheres are probably rotating very quickly. Thus, the ob- 
served photometric colors are affected by this rotation, and so 
are the derived atmospheric parameters. The photometric ages, 
determined from these atmospheric parameters and stellar evo- 
lutionary models (Asiain et al., 1997), are systematically over- 
estimated because of this effect. In order to evaluate this ef- 
fect, we have corrected the observed photometric colors for ro- 
tation by considering a constant angle of inclination between 
the line of sight and the rotation angle (i — 60°), and a con- 
stant atmospheric angular velocity (lo = 0.6 uJc, where ujc is 
the critical angular velocity). In this way we obtain a mean age 
« l.l-10'^yr,and even smaller values for cj > 0.6 uJc- More de- 
tails on the correction for atmospheric rotation can be found in 



Figueras & Blasi (1998). Second, before applying our method 
to detect moving groups (Paper I) we eliminated all stars with 
relative errors in age bigger than 100 %, which slightly bias the 
age of young groups to larger values. 
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Table 2. Number of members (A/"), spatial velocity components 
and ages (along with their dispersions) of the Sco-Cen associ- 
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Fig. 5. Evolution of the observed dispersion {dot-dashed line), 
mean propagated error (dashed line) and intrinsic dispersion 
(solid darker line) in p (distance to the Sun), determined from 
the most spatially concentrated stars in the B 1 group 



In a recent study de Zeeuw et al. (1999) carried out a cen- 
sus of nearby OB associations using HIPPARCOS astrometric 
measurements and a procedure that combines both a conver- 
gent point method and a method that uses parallaxes in addi- 
tion to positions and proper motions. Their sample of neigh- 
bouring stars is much larger than ours, since radial velocities 
and Stromgren photometry are needed for our analysis. They 
found 521 stars in the Sco-Cen association, of which only 65 
were in our sample. From these stars we have determined the 
mean velocity components of each association, which are in 
close agreement to those of Bl (Table However, their dis- 
persions and ages are quite a bit larger than expected for a typ- 
ical association. A deeper analysis revealed to us the presence 
of some stars in de Zeeuw et al.'s (1999) associations whose 



peculiar velocities are responsible for their large velocity dis- 
persion. Consistently, the ages of these stars are also peculiar 
(r « 1 — 5 • 10* yr). Removing these stars and a few others 
with anomalously large ages^ all of which probably belong to 
the field, we find a more reasonable kinematic properties and 
ages for these associations (Table 

4.2. B2 group 

Even though B2 is quite a young group, the propagation of 
age uncertainties over time prevents us from determining the 
group's kinematic age by means of the procedure developed 
for the Bl group. Instead, we have determined the trajectories 
of each star in B2, and that of the MG's nucleus itself. The spa- 
tial concentration of these stars is propagated back in time by 
counting the number of stars found in 300 pc around the MG 
nucleus (Fig. |^). We do not find any maximum in spatial con- 
centration during the last 10* yr. To better understand the struc- 



' Eliminated stars: HIP79599 and HIP80126 in US; HIP70690, 
HIP73147, HIP80208 and HIP83693 in UCL; HIP60379, HIP62703 
and HIP64933 in LCC 



ations found by de Zeeuw et al. (1999), determined from those 
members present in our sample 



Assoc. 
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U 


V 


w 


r 






[kms^M 


[kms~^] 


[kms-^] 


[10* yr] 


US 


14 


-3.5(4.3) 


-15.7(2.1) 


-6.4(1.8) 


0.4(0.5) 


UCL 


32 


-3.7(7.0) 


-19.7(4.1) 


-4.9(2.1) 


0.7(1.0) 


LCC 


19 


-6.8(4.7) 


-18.5(6.5) 


-6.4(1.7) 


0.9(1.2) 


TOTAL 


65 


-4.6(6.1) 


-18.5(4.9) 


-5.7(2.1) 


0.7(1.0) 



Table 3. Same as Table ^ now removing a few stars with pe- 
culiar velocity components and ages 



Assoc. 


u 


U 


V 


W 


r 






[kms^^] 


[kms~^] 


[kms'M 


[10* yr] 


US 


12 


-3.8(4.0) 


-15.5(2.1) 


-6.6(1.8) 


0.2(0.1) 


UCL 


28 


-1.7(4.4) 


-20.4(3.6) 


-4.4(1.5) 


0.4(0.3) 


LCC 


16 


-5.1(2.7) 


-20.6(4.0) 


-6.0(1.6) 


0.4(0.3) 


TOTAL 


56 


-3.1(4.2) 


-19.4(4.0) 


-5.3(1.9) 


0.4(0.3) 



ture of this group we plot in Fig. ^ the position of these stars and 
their velocity components on the galactic and meridional plane, 
referred to the LSR and corrected for galactic differential rota- 
tion. We observe that it is actually composed of several spatial 
stellar clumps, each of which has its own kinematic behaviour, 
although certain degree of mixture is also evident. In particular, 
one of these clumps (most of the filled circles in Fig. ^ per- 
fectly overlaps with the Sco-Cen association mentioned above 
(20 of these stars are classified as members of this association 
by de Zeeuw et al. 1999). This clumps's velocity components, 
and kinematic age, are also very similar to those found for the 
Bl group. Once again, the photometric age of these stars is 
probably overestimated because of the high rotation velocity 
of their atmospheres^ A second clump is coincident with the 
Cassiopeia-Taurus (Cas-Tau) association in both position and 
velocity spaces (most of the empty circles in Fig. 0). There are 
two dominant streams among these stars. In the heliocentric 
system, corrected for galactic differential rotation, their veloci- 
ties are (U,V)~(-10,-18) km s^^ and (-13,-23) km s~^ respec- 
tively (they share a common W component ^ —6 km s~^). 
Their ages are 5 — 6 • 10^ yr These values are in good agree- 
ment with those found by de Zeeuw et al. (1999) for Cas-Tau 
association, i.e. (U,V,W) = (-13.24,-19.69,-6.38) km s^^ with 
an age of 5 • 10^ yr (actually, about half of them are classified as 
Cas-Tau members by these authors). Finally, we observe a less 
concentrated group of stars at ^ w 0.05 kpc and 77 k, 0.2 kpc, 
or I w 90° (diamonds in Fig. The kinematic properties and 



* Our MG finding algorithm (Paper I) was able to detect 47 of the 56 
Sco-Cen association members in our sample. However, the classifica- 
tion procedure we applied could not group them into a single MG. The 
age of these very young stars, affected by the important uncertainties 
already mentioned, is probably causing this misclassification 
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Fig. 6. Number of MG stars placed around their nucleus (at a 
distance < 300 pc) as a function of time, for MGs B2, B3 and 
B4 



age of this clump are very similar to those of the Cas-Tau asso- 
ciation, whereas its position quite well agrees with that of the 
recently discovered Cepheus OB6 association (Hoogerwerf et 
al., 1997). 

Thus, B2 seems to be the superposition of several OB as- 
sociations from the Gould Belt, which are mixing with each 
other in the process of disintegration. These stars were classi- 
fied as belonging to the same MG in Paper I since they roughly 
share the same kinematics and age, a consequence of their be- 
ing formed from the same material. 

4.3. Older Pleiades subgroups 

Groups B3 and B4 in Table |l] are considerably older than Bl 
and B2, and therefore only few details on the conditions in 
which they were formed can be recovered. In particular, the 
uncertainties in 77' due to observational error increase almost 
monotonically with time. Though this increase can be closely 
approximated by Eq. || for an axysimmetric galactic potential, 
when considering the terms 4>sp and $b this approximation 
breaks (it only works during the first « 10* yr). To estimate the 
effect of typical errors in both the position (« 10-15 pc) and 
velocity (« 2-3 km s^^) components of moving groups on if 
as a function of time, we have simulated a stellar group whose 
dispersions in the phase space equals those errors, then deter- 
mined their orbits. We obtain a dispersion in t]' due to these 
errors of w 500 pc after ^ 3 • 10* yr Moreover, we cannot 
know the way in which individual orbits have been perturbed 
due to the disc heating effect, producing additional spatial and 
velocity dispersions. 

As mentioned above, these older groups seem to have been 
born in the vicinity of the spiral arms. The position at birth 
of groups B3 and B4 are especially interesting: on the one 
hand, the trajectories followed by these groups show a max- 
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Fig. 7. Position of B2 stars at present and their velocities re- 
ferred to the LSR and corrected for galactic differential rota- 
tion. Different symbols are used for three galactic longitude 
ranges. The size of the symbols is inversely proportional to the 
age of the stars they represent 



imum galactocentric distance at the moment they were born, 
which corresponds to a minimum kinematic energy (Figs. |^ 
and ||); on the other hand, by determining the trajectories of 
the B3 members we observe two focusing points close to the 
B3 and B4 birthplaces (Fig. 0). An identical result is ob- 



tained when using B4 stars. Following Yuan (1977), these 
points could be interpreted as the birthplace of B3 and B4. 
For spatially concentrated groups of stars with small veloc- 
ity dispersions epicycle theory predicts (Sect, pb the focus- 

27r 

ing phenomenon will be produced every At = — . Since 

K 

the Pleiades groups oscillate around a guiding center placed 
at i? « 8.0-8.2 kpc, the corresponding epicycle frequency is 
K « 38-40 km s^^ kpc (determined from <i>AS, Sect. 3.1), 
with At « 1.5-1.6 -10* yr This period At is compatible with 
the ages of B3 and B4. 

By following the same procedure used for the B2 group we 
study the spatial concentration of the older groups (Fig. For 
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Fig. 8. Stellar trajectories of the B3 members during the last 
3.510* yr The position of the B3 nucleus at t = 0, -1.5 and 
-3 10* yr (in units of 10* yr) is indicated with an open circle. 
The reference system is the same as in Fig. ^ 

B4 there is a maximum concentration at t w — 1.5 • 10*, which 
corresponds to this MG's age, whereas we observe two peaks 
in concentration for B3 at t w —1.5 and —3 • 10* respectively, 
the first one corresponding to the last focusing event, and the 
second one corresponding to the average age of its members. 

5. Evolution of a Stellar Complex 



As defined by Efremov ( |1988[ ), Stellar Complexes (SC) are 
"groupings of stars hundreds of parsecs in size and up to 
10* yr in age, apparently uniting stars born in the same gaseous 
complex". Associations and open clusters are the brighter and 
denser parts of these huge gaseous complex es. In other words, 
and according to Elmegreen & Elmegreen ( 1983 ), the funda- 
mental unit of star formation is an Hl-supercloud of ~ 10^ Mq 
which, during fragmentation into ^ 10^-10^ Mq giant molec- 
ular clouds, gives birth to a SC. From these giant molecular 
clouds clusters and associations are formed. Examples of SCs 
are the Gould Belt in our Galaxy, 30 Doradus in the LMC, and 
NGC 206 in Andromeda. 

In this section we explore the possibility that such objects 
are the progenitors of MGs. Since both galactic differential ro- 
tation and disc heating effect tend to disrupt any unbound sys- 
tem of stars on a short timescale, the large number of stars born 
inside a SC - and therefore roughly sharing the same kinemat- 
ics and ages - may ensure a high spatial concentration for a 
longer time, especially at the focusing points. 

In order to analyse the possible link between SCs and MGs, 
we have generated a SC as a set of different unbound systems 
born at different epochs. Its main characteristics are taken from 
the literature and are described in what follows. Although there 
are no preferences on the galactic plane region where SCs are 
formed (they can even be found in the interarm region, as sug- 



gested by the observations of OB Associations in Andromeda, 
e.g. Magnier et al. 1993), we select the position and velocity 
of the B4 nucleus at birth (t — —1.5 • 10* yr), since in this 
way we can compare the results here obtained with observa- 
tions on this real group. According to several authors' estima- 



tions (Efremov, 1988; Elmegreen, 1985; Efremov & Chernin, 



1994), SCs are several hundreds of parsecs in size. The orig- 
inal shape of our simulated complex has been designed as an 
ellipsoid whose equatorial plane is 250 pc in radius and lies on 
the galactic plane, and its vertical semi-axis is 70 pc. Six as- 
sociations are born inside it at different epochs, randomly dis- 
tributed inside this volume. Since the disp ersion in age ca n be 
large when considering the SC as a whole ( Efremov, 1988 ), we 
impose the condition that one associations borns every 10^ yr. 
Hence, the first one is born at i = — 1.5 • 10* yr, the second one 
at —1.4 • 10* yr, and so on. The mean velocities of these associ- 
ations at birth follow a gaussian distribution around the B4 nu- 
cleus velocity, with an (isotropic) dispersion similar to the ve- 
locity dispersion inside a molecular cloud, i.e. w 5/\/3 km s^^ 
in each component (Stark & Brand, 1989). Each of these asso- 
ciations is defined as follows: they contain 500 stars, which are 
randomly distributed inside a sphere of 15 pc in diameter, an 
intermediate value between the smallest associations (e.g. the 
Trapezium cluster) and the largest ones (Blaauw, 1991; Brown, 
1996); and their internal velocity dispersion is 2 km s^^ in each 
velocity component, as observed in most of the closer OB As- 



sociations (Blaauw, 1991) and in molecular clouds (Scoville, 
1990); their internal dispersion in ages is 10 '' yr (Efremov, 
1988). We also assume that all the simulated stars survive up to 
the time (t) when we study the system. 

The orbit of each simulated star from the moment it was 
born up to the present (t = 0) is determined considering the 
galactic gravitational potential ($as + ^'Sp + '^'b). without 
taking into account the disc heating effect. Their current spa- 
tial distribution is shown in Fig. |[ where different symbols are 
used for the stellar "groups" coming from each simulated as- 
sociation. According to our simulation, these groups are still 
concentrated in space, although their distribution is, at present, 
far from being isotropic due to the differential rotation. Each 
of them occupies an ellipsoidal volume « 150-250 pc wide, 
500 pc long, and 70 pc height, and are perfectly separated in 
space (four of them can be observed at different positions inside 
a 300 pc radius sphere around the Sun). The dispersion in each 
velocity component of any of the groups is still « 2 km s^^, as 



expected from Eqs. |A.4 Although the mean velocity of these 
stars is very similar to that of B4, the dispersions in both ve- 
locity and space are too small when compared to dispersions 
in detected MGs that ai-e a few 10* yr old (Table |l|). Thus, the 
spiral arms and central bar perturbations to the galactic grav- 
itational potential cannot account for the observational MGs' 
velocity dispersions. 

Results drastically change when considering the effect of 
constant heating on each individual star's trajectory, as de- 



scribed in Sect. 3.2. To determine these trajectories we use only 



the axisymmetric part of the galactic potential, since the asym- 
metric perturbations to this potential are already included in 
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Fig. 9. Position on the galactic (a) and meridional (b) planes of the simulated stellar complex at t = (see text) when no disc 
heating effect is considered. Three square regions (600 pc/side) on the galactic plane have been enlarged to better distinguish 
the details. A different symbol is used for each group of stars coming from the same association. GC: galactic center; NP: north 
galactic pole; GR: galactic rotation 



the treatment of the heating. The distribution of the simulated 
stars at present is shown in Fig. |lO| The members of the dif- 
ferent groups are now completely mixed with each other The 
total SC forms a long structure ~ 3 kpc long and ^ 1 kpc wide. 
Since in this simulation the Sun has been considered to be in the 
very center of the SC at < = 0, the members of this complex 
can be found up to some 500 pc from us in the radial direc- 
tion, and much further in the tangential direction. By enlarging 
a square region (600 pc/side) around the Sun we observe in 
Fig. |lo|that, although the members of two groups dominate in 
the solar neighbourhood, all the other evolved associations are 
also represented. Again, the mean velocity components of B4 
are recovered when considering only the (636) stars closer than 
300 pc from the Sun. However, as expected, the dispersions in 
velocities are now higher than in the former simulation - the 
total velocity dispersion is w 8-10 km s^^, depending on the 
region of the evolved SC. Hence, a constant diffusion of the 
stellar velocities can perfectly account for B4 velocity disper- 
sions. 

In order to study the properties of younger MGs, we anal- 
yse the general properties of the complex at its earlier stages, 
only 5 • 10^ yr after the first association was born. At that mo- 
ment (t = — 1 • 10^ yr) some stars in the youngest simulated 
groups are not yet born, while some others are as young as B 1 



and B2 members. The associations are still very concentrated 
in the phase space, though some merging can be observed, as 
happens with B2 group. Velocity dispersions inside each asso- 
ciation clearly depend upon their ages, as expected. Thus, in the 
U-component the velocity dispersion varies from w 2.7 km s~^ 
for youngest groups to « 4. 1 km s^^ for the oldest (in the other 
components the dispersions are smaller). 

Finally, it is interesting to take a look to the general prop- 
erties of a simulated SC at a larger t. The mean properties of 
the simulated SC are now taken from the B3 nucleus at birth 
(t = —3 • 10^ yr). At present (t = yr), assuming a constant 
heating of ah = 1-45 km s^^ every At = 10'' yr, the 3000 
simulated stars occupy a huge curved region, about 2 kpc wide 
and almost 10 kpc long. In the most dense parts of this region, 
~ 200 stars can be found in a 300 pc radius sphere, whose total 
velocity dispersions are very high (12-14 km s^^). Stars with 
such large dispersions would be completely merged with field 
stars, so they could not be detected as MG members. A much 
smaller diffusion coefficient is needed to recover the B3's ve- 
locity dispersions (ah ~ 0.7 km s~^). 

Thus, the disgregation of SCs by means of a constant heat- 
ing mechanism is not able to account for the whole set of 
Pleiades MG substructures. A diffusion coefficient that de- 
pends on the stellar peculiar velocity and/or on the time (e.g. 
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Fig. 10. Same as Fig. ^but considering now the effect of the disc heating on the trajectory of the stars 



episodic diffusion) probably might explain the observed prop- 
erties of these groups. Moreover, the presence of some open 
clusters in the original SC could keep the stars together in phase 
space during longer periods. Some young open clusters, i.e. 
IC 2391, a Persei, Pleiades, etc, share roughly the same kine- 
matics as the groups in Table [l|, which favours this hypothesis. 

6. Conclusions 

In this paper we studied the disruption of unbound systems of 
stars as a mechanism to understand the origin and evolution of 
moving groups. The epicycle theory was used to find an ana- 
lytic expression for the time dependence of dispersion in both 
the stellar position and velocity coordinates. We have obtained 
in this way a simple expression for the secular increase of the 
dispersion in the azimuthal galactic coordinate over time as a 
function of the initial conditions. This increase in dispersion is, 
in fact, a direct consequence of the galactic differential rotation. 

In order to overtake the constraints of the epicycle theory, 
we concentrated our analysis on the determination of the stellar 
trajectories using numerical integration of the equations of mo- 
tion, which provided us with an independent and more accurate 
estimation of the evolution of unbound systems. To perform 
this integration we used an analytic and axisymmetric galactic 
potential, along with spiral arm and central bar perturbations. 
This latter procedure allowed us to include random perturba- 
tions that mimic the disc heating effect on stellar trajectories. 



The trajectories followed by the members of the Pleiades 
moving group substructures found in Paper I were used to com- 
pare the kinematic and photometric ages of these structures, 
and to establish their position at birth. The youngest group, Bl 
(the Sco-Cen association), was found to be most spatially con- 
centrated some 9 • 10® yr ago, a value considerably smaller than 
its photometric age. This is probably due to the effect of the 
high stellar atmospheric rotation of early type stars on the ob- 
served photometric colors. Groups B4 and B3 were born at a 
time from the present equal to one and two times the epicy- 
cle period respectively, which means that they are spatially fo- 
cused at present (probably that is the reason why we can ob- 
serve them). At birth, these two groups were found to be close 
to the spiral arm structure. Concerning B2, a detailed analysis 
revealed to us that it is actually composed of several associa- 
tions which are disintegrating at present. As well, its averaged 
photometric age is probably overstimated, as in the case of B 1 . 

We considered the evolution of a stellar complex under the 
influence of the galactic gravitational potential as a mechanism 
to account for the main physical properties of moving groups. 
The high velocity dispersions of some of the Pleiades moving 
group substructures detected among B and A type stars could 
be recovered when the effect of the disc heating on the indi- 
vidual stellar trajectories was considered. At the same time, 
the disc heating can account for the mixing of the stellar com- 
plex associations, although they are still clearly separated in 
phase space during the first tens of million years in the com- 
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plexes' life (as observed for the groups Bl and B2). After only 
1.5 • 10^ yr from the birth of the stellar complex, the complex 
occupies an ellipsoidal area of ^ 3 x 1.5 kpc^, with its longest 
axis oriented in the direction of galactic rotation. If the Sun 
were presently located in the very center of this disrupted stel- 
lar complex, we would be able to find the complex's compo- 
nents up to very large distances on the galactic plane. 

Thus, a constant heating mechanism (compatible with the 
observational heating law) acting on the stars of a stellar com- 
plex can explain the velocity dispersions obtained for those 
Pleiades moving group substructures that are younger than 
~ 1.5 • 10^. The properties of the older Pleiades substructures 
could probably be recovered by considering a diffusion coef- 
ficient depending on the velocity of the stars, and maybe on 
the time, and/or the inclusion of open clusters in the simulated 
stellar complex. 
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Appendix A: Propagation of dispersions with time in the 
epicycle theory 

Let us consider the vector x = ,rj' X' j)' X') containing 
the position and velocity of a star at a given time t. We 
are interested in determining the propagation of the obser- 
vational dispersions over time, in the epicycle approximation 
frame (Eqs. [l| and their derivatives). If we define the vector 
^° = ("^o'^o'Co'Co'^o'Co) which contains the initial position and 
velocity of our star, then the initial dispersions can be propa- 
gated with time as: 



dx° 



(A.l) 



where ax° is the observed dispersion in the variable Xi, 
and <Jxi{t) the dispersion at t. The dependence of variables 
X on the initial values x° is given by Equations |[ If the 
coefficients of epicycle equations are included in a vector 



i^UM',^^^) then the matrix A = (a)) = 1^ 



in 



Equation A. 1 can be determined as the matrix product A = BT, 
where the matrix B = {h)) = and T = (/]) = 



dcj 
dx° 



1, 



, 6). By partially differentiating Eqs. |] 



we can determine the matrix B (Table A. 1 ), whereas matrix 
is determined by just differentiating Eqs. 2| (Table A. 2). Now, 
from matrices B and T we can determine the matrix A, given 
in Table A. 3 From matrix A, and by means of the relation- 
ship A. 1 , the propagated dispersions can be expressed as: 
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1 ^ 


sin(Kf) 


al + 




K 



4 + 



2B 



2tJoA f sm{Kt) 
B ■ ^ 

2u)„ 



1 2 



(cos(Kt) — 1) 



— (At — — sm(Kt) 
B V K 



4 

So 

n 2 



2 2/ ,\ 2 I 

(T^, = COS {vtjCjQ,^ + 

' Ak 



B 

K 

2B 

2Acjc 
B 



sin(Kt) 



sin(z^i) 



fj?/ + cos'^ (Kt)aj, 

so 'to 



4. 



sin(Ki) 
■ (1 — cos{Kt)) 



1 2 



2Wn 



1 2 



■ sin(Ki) 



4 

So 



— (A — Wo cos(Kt)) 
B 



(A.2) 



From these equations we observe that all the dispersions 
in variables x oscillate around a constant value, except cr^,, 
whose averaged value increases with time. For small values of 



t Eqs. A.2 can be approximated by: 



^2 , .2 2 
Gi, 

So 

So 
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